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Abstract 

A model for a homogeneous and isotropic spatially flat Universe, composed of baryons, radiation, 
neutrinos, dark matter and dark energy is analyzed. We infer that dark energy (considered to behave 
as a scalar field) interacts with dark matter (either by the Wetterich model [1], or by the Anderson 
^^^^ and Carroll model [2]) and with neutrinos by a model proposed by Brookfield et al. [3]. The latter 

I I is understood to have a mass-varying behavior. We show that for a very-softly varying field, both 

{2JQ' interacting models for dark matter give the same results. The models reproduce the expected red-shift 

performances of the present behavior of the Universe. 

PACS: 98.80.-k, 98.80.Jk 

O ' 1 Introduction 

rn '. 

^T) [ Astronomical data observation, such as measurements of type-IA supernovae [4, 5, 6, 7] indicates that our 

O ' Universe is ongoing an accelerated expansion period, suggesting the existence of a different component, 

I ■ with the strange propriety of negative pressure, whose energy density is now predominant over others. This 

^* i entity was named dark energy, and its nature still remains a mystery. Besides, the measurement of rotation 

• ^H ' curves of spiral galaxies [8] is not in agreement with the luminous matter within it, proposing a new kind of 

rS ' matter, which does not radiate, and interacts with common components only gravitationahy, named dark 

d . matter. Since those discoveries, researches have been made concerning these constituents behavior. Dark 

energy may be treated as the cosmological constant [9], but may also have dynamical proprieties, inciting 

studies with several kinds of scalar fields. Some earlier studies considered the interaction of a scalar 

field with any kind of particles [10], suggesting the further treatment of particular cases. Furthermore, 

some puzzling questions arise from observation, such as the fine-tuning and the coincidence problem 

[11, 12, 13, 14, 15], giving rise to a whole new investigation of the components of our Universe, in which 

dark energy shall interact with both dark matter and neutrino [1, 2, 3, 16, 17, 18, 19, 20, 21]. The 

latter has the great advantage of experimental study, even though neutrino interaction is quite weak, the 

abundance of this particle allows not only detection, but opens road to investigation of its oscillating mass 

evolution, permitting one to develop more precise theoretical models, defining bound limits for its mass. 

In this paper, we work with the spatially flat Friedmann-Robertson- Walker metric, and analyze the energy 

densities evolution under two different dark energy-dark matter coupling, the Wetterich interaction [1, 22] 

and the Anderson and Carrol interaction model [2, 23], which we find out to be equivalent. Neutrinos are 

considered to present a mass-varying behavior and to be coupled to dark energy by a model proposed by 

Brookfield et al. [3]. Natural units were rescaled (SvrG = c = 1). 
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2 Field Equations 

We consider a Universe composed by the following fields: baryons (6), radiation (r), dark energy (99), 
dark matter (dm) and neutrinos (v). The baryons and radiation are regarded as non-interaction fields, 
the dark energy is interpreted to behave as a scalar field which interacts with the dark matter and the 
neutrinos are treated as a mass- varying component due to a coupling to the dark energy field. 

In a conformal time the Friedmann-Robertson- Walker metric for an isotropic, homogeneous and spa- 
tially fiat Universe is given by 



ds^ = airy dr' - 6ijdx'dx^ , (1) 

where a{T) represents the cosmic scale factor. 

The Friedmann and acceleration equations which follows from Einstein field equations in conformal 
time read 

a/ 6 a b 

respectively. The dots denote the differentiation with respect to the conformal time, while p = Pb -\- Pr ~\~ 
Pdm + Pu + Pip and p = Pb + Pr + Pdm + Pu + Ptp are the total energy density and pressure of the sources, 
respectively. 

From the conservation of the energy-momentum tensor it follows the evolution equation for the total 
energy density: 

p + 3-(p + p)=0. (3) 

a 

The baryons are treated like dust {pb = 0) and the barotropic equation of state for the radiation field 
is Pr = Pr/'i- Once both of them have no interaction with other components, the evolution equations for 
their energy densities are: 

Pb + 5-pb = 0, Pr + ^-Pr = 0, (4) 

a a 

respectively. 

Since the dark energy is modeled as a scalar field if its energy density and pressure are given by 

1 1 

P'P = ^^^ + y(.y^), Pv = i^v>^ -y{^), (5) 

where V{ip) denotes the potential of the scalar field. 

In this work we shall analyze two types of interactions between dark matter and dark energy. The 
first one is the Wetterich model [1] which for a pressureless {pdm = 0) interacting dark matter has the 
following evolution equation 

Pdm + ^-Pdm = -PPdm'P- (6) 

a 
The other refers to the Anderson and Carroll model [2] and in this case the evolution equation for the 
dark matter reads 

Pdm + "i-Pdm = (3pdm — - (7) 

a (p 

In both equations above, f3 stands for the coupling constant between dark matter and dark energy, which 

may be interpreted as the fundamental scalar charge q of the particles [24] . 

The neutrinos are understood to be massless and relativistic particles in the past, but with the coupling 

to the dark energy, they have acquired mass and became non-relativistic, having an oscillating mass 

behavior at low red-shifts. The evolution equation for its energy density according to Brookfield et al. [3] 

is 

d , din nil, , \ , \ 

Pu + 'i-{Pu+Pu) = — 3 (p{pu-?>Pu), (8) 

a dip 

where niu is the neutrino mass which is considered a function of the scalar field ip. Here, we follow [3] 
and write the neutrino mass as rrii, = Mqc""^, where Mq and a denote coupling constants. 



From the evolution equation for the total energy density (3) and by taking into account (4), (8) and 
(6) or (7), we may obtain two evolution equations for the energy density of the scalar field, one referring 
to the Wetterich model and the other to the Anderson and Carrol model. Furthermore, from the resulting 
equations it is possible to obtain evolution equations for the scalar field by considering the representations 
for its energy density (5)i and for its pressure (5)2. These equations reads 



„ . a 2 

(p + Zip — h a 
a 



ip + 2ip) — h o 
a 



dy d In nil, 
dip dip 
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0, 



(9) 



(10) 



where the first one refers to the Wetterich model, while the second one to the Anderson and Carrol model. 
From now on we shall use the red-shift z as variable instead of the conformal time r. Since the red-shift 
is represented in terms of the cosmic scale factor by a = 1/(1 + z), we have the following relationships: 



d_ 



fp d 
3 dz^ 



dT^ 



p + p 1 d p d"^ 

2 l + zd^^3dz^' 



(11) 



The differential equations for baryons (4)i and radiation (4)2 can be integrated and this yields the 
well known relationships: 

Pk{z)=p^,{l + zf, pr{z) = p',{l + Z)\ 



where the index stands for the present value of the variable. 

For the energy density of the dark matter, we obtain through the integration of (9) and (10) 
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Al + zYe 
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(12) 

(13) 
(14) 



for the Wetterich and Anderson and Carroll models, respectively. 

In terms of the red-shift the differential equation for the energy density of the neutrinos (8) becomes: 



Pu 



Pu+Pu , dlnrriu, 
3-r- \ ; [Pu - 3Pv 



l + z 



dz 



(15) 



where the prime denotes differentiation with respect to the red-shift (z). 

Furthermore, the differential equations for the scalar field in terms of the red-shift for the Wetterich 
and Anderson and Carroll models - which follow from (9) and (10) - read 
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(16) 
(17) 



respectively. 

Due to the fact that the energy densities of the baryons, radiation and dark matter are determined by 
the expressions (12) - (14), the equations (15) - (17) compose our system of coupled diff^erential equations 
for pi, and 99 we are willing to solve. 

Once the neutrino mass was supposed to have the form rui, = Mqc""^, we have to specify the potential 
of the scalar field V{ip) and the equation of state of the neutrinos Pi, = Pu{pu) in order to solve the system 
of coupled differential equations. 

Since authors began to study neutrinos coupled to dark energy, many kinds of pressure have been 
attributed to it. Some of them suggest that the pressure would be proportional to their density through 
a barotropic equation of state p^ = w^Pu^ where Wu is a constant. Note that if Wy = 1/3, the coupling 
vanishes making neutrinos non-interacting particles. Recently, the neutrino varying-mass became a strong 



aim of study among researchers and this emerges the need of investigating varying-pressures, so that we 
could better explain the coupling mechanism. Here, we introduce two possible equations of state: 

in accordance with what we postulated before. We should notice that as the red-shift tends to high values, 
these equations of state turn out to be equivalent to the case of radiation [jp^ = Pi//3) during the period of 
non-coupling with the dark energy, while neutrinos were still relativistic. In the recent era the neutrinos 
are pressureless {py = 0) and they are treated as dust particles. Note that the second equation of state 
tends more rapidly to {py = Pu/'i) for large values of the red-shift z. Although we have proposed the 
two above equations of state for the description of the past and present behavior of the Universe, one 
may observe that the extrapolation to negative values of the redshift will lead to negative values of the 
neutrino pressure field and even the first equation will diverge at z = —1. Hence, the two equations of 
state proposed here cannot be extrapolated to negative values of the redshift. 
The potential adopted in this work is the "Mexican hat" potential: 

V{ip) = \[p{zf-^t\ -\^l (19) 

where A and (pQ denote two constants. 

3 Solutions 

The numerical solutions of the system of equations (15) and (16) and of (15) and (17) - according to the 
Wetterich and Anderson and Carrol models - were found by considering the equations of state (18) and 
the specification of the coupling and potential parameters together with the initial conditions. 

First we introduced the density parameters $7^ = pi/ p and their values at z = 0: ^^(O) = 0.05, 
Qr(0) = 5 X 10~^ and 0(^^(0) = 0.23. Next the initial conditions for 99, ip' and py were specified at 
z = 3100, namely, (/3(3100) = ipo = 2.69 and (/9'(3100) w and Oi.(3100) sa 0.68 0^(3100) (see e.g.[25]). 
Furthermore, as in any phenomenological theory, the investigation of the interactions in the dark sector 
requires the knowledge of some values for the coupling constants, in order to reproduce the observed 
behavior of the present Universe. Here we have the freedom to choose the potential constant and the dark 
matter coupling parameter. The potential constant was chosen to be A = —0.015, and the dark matter 
coupling parameter: /3 = —0.01 in the Wetterich model and j3 = 0.03 for the Anderson and Carroll 
interaction. 

The scalar field, for both models, varies pretty slowly during the expansion of the Universe (see Fig. 
1), which means that the energy density and the pressure of the dark energy have always been dominated 
by its potential and, except for its tiny variation, it behaves like a cosmological constant. For this reason, 
the coupling between neutrino and dark energy must be strong. Indeed, such couplings are very unusual 
and not easy-to-swallow, once free cosmological parameters are often found in a range [—1,1], suggesting 
a re-scale of the scalar field: (p = IQip + K where K = —23.83 is a constant chosen to keep the initial 
conditions. The determination of the coupling parameter a was found through the variation of its value 
and by constraining that the neutrinos mass shall present an oscillating evolution, as well as assuring 
that the dark energy density parameter must not grow at large red-shifts. Once we did this, our coupling 
parameter to neutrinos is a = 10, still a high value, but acceptable considering the field behavior, and 
the investigation of new forces. 

The scalar field varying very softly during all times, exactly the same way for both models, and not 
presenting singularities, resulted in the pretty same curves for all four models (considering now the two 
barotropic equations for neutrinos). Nevertheless, we notice from Figure 1 an oscillating behavior in 
the scalar field, which increases its value as the redshift approaches zero. Thus we may already expect 
a similar behavior for neutrinos mass. We can see in Figures 2 and 3, the evolutions of the density 
parameters, which occurred as expected: at high red-shifts, dark energy is neglected, the radiation and 
baryons fall with a~^ and a~^ respectively. We note from Fig. 3 that the density parameter of the dark 
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Figure 1: Scalar field vs. red-shift in logarithmic scale for a = 10 and Mq = 10" 
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Figure 2: Density parameters vs. red-shift. 



matter presents a very similar behavior as the baryons. This fact can be understood as follows: besides 
the energy transfer of the dark components, dark matter also falls with a~^, once the amount of dark 
energy is so small for the period shown in this figure. Furthermore, as can be observed from Fig. 2 the 
evolution of the dark energy density does not show up in Fig. 3, because it decreases very rapidly with 
the red-shift. In this model the transition between two important epochs of our Universe, the radiation- 
matter equality {i^dm + ^b = ^r + ^u) takes place at z ~ 3400. The neutrino showed a moderate energy 
density in past, decreasing very softly until recent times so that i7jy(0) ~ 1.4 x 10~^. For small red-shift 
values, the dark energy becomes dominant very fast, reaching 0(^(0) ~ 0.71 and making the Universe 
to expand acceleratedly, as we can see in Fig. 4, where we plotted the deceleration parameter, given 
by: Q = 2 "I" 2n' ^^^ z = the deceleration parameter has the value q{0) ~ —0.53, while the transition 
decelerated-accelerated occurs at zt ~ 0.74, both values are of the same order of magnitude of the ones 
given in the literature: q{0) = -0.46 ± 0.13 (see [26]) and zt = 0.74 ± 0.18 (see [27]). 

The mass evolution is plotted in Fig. 5, the constant Mq in the expression rrii, = Moe"*^ plays an 
essential role as a free parameter, and here we chose it to be Mq = 10~^^. It is clear that it does not 
change the shape of the curve (the neutrinos mass behavior does not depend on it) but instead, changing 
its value translates the curve upwards or downwards. We see from this figure that for recent times, the 
neutrino mass would be about 0.31ey. If we go back in time, as the red-shift increases, so does the 
neutrino mass, which evolves until z ^ 3000 reaching ttt-j^ = 0.21ey and then begins to decrease again 
until zero at large red-shifts. 
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Figure 3: Evolution of density parameters in large scale. 




Figure 4: Deceleration parameter vs. red-shift. 




Figure 5: Evolution of neutrino mass vs. red-shift in logarithmic scale. 



4 Conclusions 

Although we chose to work with the "Mexican hat" potential, other kinds like the Yukawa potential have 
been investigated, giving the same soft variation of the field. For this reason, for any red-shift value, 
93' assumes no predominant position, making the field pressure and energy density dominated by its 
potential. Furthermore, the almost constant behavior of the field, makes both interacting models differ 
only by a constant, which is contoured by attributing correct values to the coupling parameters, resulting 
in equivalent models. As mentioned before, it suggests a very strong dark energy-neutrino coupling, 
if we are really attained to neutrino mass-varying property, leading to the assumption that the deep 
investigation of the latter will improve our understanding of the first. 
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